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Introduction: The dust cycle is critically 

important for the current climate of Mars. The 
radiative effects of dust impact the thermal and 
dynamical state of the atmosphere (Gierasch and 
Goody, 1968; Haherle et al., 1982; Zurek et al., 
1992). Although dust is present in the Martian 
atmosphere throughout the year, the level of 
dustiness varies with season. The atmosphere is 
generally the dustiest during northern fall and winter 
and the least dusty during northern spring and 
summer (Smith, 2004). Dust particles are lifted into 
the atmosphere hy dust storms that range in size 
from meters to thousands of kilometers across 
(Cantor et al., 2001). During some years, regional 
storms combine to produce hemispheric or planet 
encircling dust clouds that obscure the surface and 
raise atmospheric temperatures by as much as 40 K 
(Smith et al., 2002). Key recent observations of the 
vertical distribution of dust indicate that elevated 
layers of dust exist in the tropics and sub-tropics 
throughout much of the year (Heavens et al., 2011). 
These observations have brought particular focus on 
the processes that control the vertical distribution of 
dust in the Martian atmosphere. The goal of this 
work is to further our understanding of how clouds 
in particular control the vertical distribution of dust, 
particularly during N. H. spring and summer. 

The past few years have seen an increase in the 
recognition of the importance of water ice clouds on 
the state of Mars' current climate. The main surface 
source of atmospheric water is the north residual 
cap. Broadly, there are two categories of clouds: the 
aphelion cloud belt and the polar hood clouds. The 
aphelion cloud belt is composed of optically thin 
clouds that fonn above 10-15 km at low latitudes 
during northern spring and summer (Ls~50°-135°; 
Clancy et al., 1996). Polar hood clouds form near 
the edge of the advancing and receding polar cap in 
both hemispheres during local fall and winter. 
These clouds are optically thick and may or may not 
extend down to the surface (Clancy et al., 1996; 
Benson et al., 2010; Benson et al., 201 1). 

Cloud formation is the key process for the 
coupling between the dust and water cycles (Figure 
1). Dust particles likely provide the seed nuclei for 
heterogeneous nucleation of water ice clouds 
(Montmessin et al., 2002; Maattanen et al., 2007). 
As ice coats atmospheric dust grains, the newly 


formed cloud particles exhibit different physical and 
radiative characteristics. As the mass ratio of ice-to- 
dust increases, the particle’s scattering properties 
become those of ice instead of dust. The mass and 
size (and therefore bulk density) of the cloud 
particles depend on the fraction of ice versus dust. 
Because the gravitational fall velocity of particles 
depends on these quantities, cloud formation will 
either increase or decrease the suspension lifetime of 
atmospheric dust. The coupling between the dust 
and water cycles most likely affects the atmospheric 
distributions of dust, water vapor and water ice and 
thus atmospheric heating and cooling and the 
resulting circulation and dynamics. We are 
interested in how these couplings affect the dust 
cycle itself 



Figure 1 : Key processes of Mars’ dust and water cycles. 


The NASA Ames GCM: The NASA Ames 

GCM is a 3D finite -difference model of the Martian 
atmosphere that has been used extensively for 
investigations of Mars' current and past climate 
(Haberle et al., 1999; Kahre et al., 2006; Kahre and 
Haberle, 2010; Hollingsworth and Kahre, 2010). 
The model runs on a latitude/longitude horizontal 
grid and a nonnalized sigma-coordinate vertical grid. 
For this study, the nominal horizontal resolution is 
5° by 6°. Surface properties include MOLA 
topography and albedo and thermal inertia maps 
derived from Viking and MGS/TES observations. 
This version of the model includes a 2-stream 
radiative transfer scheme that accounts for gaseous 
absorption and scattering aerosols and a tracer 
transport scheme based on the Van Leer formulation. 

Routines have been incorporated into the GCM 
to account for the physics of the lifting, transport and 
sedimentation of radiatively active dust. Two 
parameterizations for dust lifting, the wind stress and 
dust devil schemes, are used simultaneously. The 
wind stress lifting scheme lifts dust when the 


momentum imparted to the surface (characterized as 
the surface wind stress, x) exceeds a critical value, 
assumed here to he x=22.5 mN m'^. The dust devil 
lifting scheme is based on the thermodynamic theory 
of dust devils developed hy Renno (1998). In this 
scheme, the lifted dust flux depends on the 
magnitude of the sensible heat exchange between the 
surface and atmosphere and the depth of the 
planetary boundary layer. The lifted dust mass is 
partitioned into a lognormal distribution with an 
effective radius of 2.5 pm and an effective variance 
of 0.5. Airborne dust interacts with solar and 
infrared radiation, provides seed nuclei for water ice 
clouds, and undergoes gravitational sedimentation as 
free dust and as cores of water ice cloud particles. 
Both dust lifting scheme are tuned with 
multiplicative “efficiency” factors to produce 
reasonable dust loadings throughout the Martian 
year (see, for example, Kahre et al, 2006 for a 
discussion about tuning). 

The simulated water cycle includes sublimation 
from the north residual cap and the microphysical 
processes of nucleation, growth, and settling of 
water ice clouds (Montmessin et al., 2002, 2004). 
The particle size distribution of cloud particles is 
represented in the model using a “momenf’ method, 
whereby it is assumed that the distribution is 
lognormal and therefore fully described by a 
spatially and temporally varying mass and number, 
and a fixed (constant in time and space) effective 
variance. This allows for the evolution of cloud and 
dust particle sizes as the result of microphysical 
processes in a computationally efficient manner. 
Water ice clouds can either be radiatively active or 
inert. 

Two simulations are presented: the first includes 
an interactive dust cycle that is fully coupled to the 
water cycles through cloud formation; the second 
includes an interactive dust cycle that does not 
include cloud formation. In both simulations, the 
dust lifting schemes are tuned identically (i.e., they 
have the same efficiency factors), which allows for 
the most straightforward comparison of the results. 

Results: Although both simulations are multi- 
annual, only results from the first half of one 
simulated year is shown. Focus is given to the first 
half of the Martian year in order to quantify the role 
that clouds in the aphelion cloud belt play in 
controlling the vertical distribution of dust. 

Seasonal Patterns of Dust and Water Ice Clouds: 
In both simulations, the dust devil lifting scheme is 
responsible for providing the majority of the 
atmospheric dust during the first half of the year (Lj 
0°-180°). This result is consistent with several 
previous investigations (Newman et al., 2002; Basu 
et al., 2004; Kahre et al., 2006). Column 9-pm dust 
opacities range from zero to a few tenths, which is in 
general agreement with MGS/TES-observed dust 


opacities during these seasons (Figure 2). Increased 
opacities occur along the retreating edge of the 
seasonal north CO 2 ice cap as the result of increased 
wind stress lifting during early N.FI. spring, 
particularly in the simulation that includes cloud 
formation. Radiatively active polar hood clouds in 
the north increase the baroclinic eddy activity along 
the growing and receding north CO 2 cap, which 
leads to increased dust lifting (Flollingsworth et al., 
2011; Kahre et al., 2011; Wilson et al., 2011). 
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Figure 2: Zonal mean column 9-pm dust opacities during 
northern spring and summer (Ls=0°-180°) for the 
simulation without clouds (top) and the simulation with 
clouds (bottom). 
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Figure 3: Zonal mean column 12-pm water iceopacities 
during northern spring and summer (Ls=0°-180°) for the 
simulation with clouds. 

In the simulation with cloud formation, the 
aphelion cloud belt forms just after the northern 
vernal equinox, dissipates just before the northern 
autumnal equinox, and spans from approximately 30 


S to 30 N. Column 12-|im ice opacities reach 
approximately 0.2 slightly to the north of the equator 
(Figure 3). The spatial extent, seasonality, and 
column opacity of the simulated aphelion cloud belt 
is generally consistent with MGS/TES observations 
(Smith, 2004). Although retrievals of column ice 
opacities over the seasonal CO 2 ice caps is difficult, 
comparing the opacities of the edges of the 
simulated polar hoods to available TES and MCS 
data suggests that the model is over-predicting the 
thickness of the polar hoods. 

Column-Integrated Dust and Water Ice Cloud 
Opacities at 1^-120°: Predicted daily mean local 

column dust IR opacities at Ls~120° range from zero 
to -0.5 in the simulations with and without cloud 
formation (Figure 4). In both simulations, column 
dust opacities maximize in the Northern Flemisphere 
over the low thermal inertia continents of Tharsis, 
Arabia, and Elysium. Simulated dust devil lifting is 
particularly active in these regions, in part due to the 
large temperature difference between the surface and 
the near-surface atmosphere (and thus the sensible 
heat flux). The most notable differences between the 
two simulations occur in the Southern Flemisphere. 
The atmosphere over the south CO 2 seasonal ice cap 
is clearer of dust in the simulation with cloud 
formation than in the simulation without, which is 
directly attributable to the thick polar clouds that 
form in the former case that effectively scavenge 
dust from the atmosphere (Figure 5). 
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denotes the edge of the seasonal CO 2 ice cap. 

Predicted daily mean local column water ice IR 
opacities at Ls~120° range from zero to >2.0 in the 
simulation with cloud formation (Figure 5). 
Optically thick clouds form near the surface over the 
south seasonal CO 2 ice cap. In general agreement 
with observations and previous modeling studies, the 
aphelion cloud belt is zonally asymmetric, with 
maxima predicted over the high topographic relief 
regions of Tharsis, Syrtis and Elysium (REFS). 
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Figure 5: 10-sol mean column 9-pm dust opacities 
Ls=120° for the simulation without clouds (top) and the 
simulation with clouds (bottom). The solid yellow line 
denotes the edge of the seasonal CO? ice cap. 

Vertical Distribution of Dust and Water Ice 
Clouds at 15 N at 1,-120°: We focus our analysis at 
15 N latitude during the afternoon (i.e., 12-3 pm 
local time) in order to understand how water ice 
clouds and their radiative effects in the aphelion 
cloud belt affect the vertical distribution of dust at 
low latitudes during N. H. summer. 
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Figure 4: 10-sol mean column 9-p.m dust opacities 
Ls=120° for the simulation without clouds (top) and the 
simulation with clouds (bottom). The solid yellow line 
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Figure 6: 10-sol mean dust (shaded) and cloud mixing 
ratios (contoured) at Ls=120° for the simulation without 
clouds (top) and the simulation with clouds (bottom). 

The predicted water ice clouds reside in the 
afternoon between 1 and 0.1 mhar, with maximum 
mass mixing ratios occurring at approximately the 
0.5 mhar pressure level (Figure 6, bottom panel). 
These simulated clouds are lower in the atmosphere 
than MRO/MCS observations suggest at this season 
(Heavens et al., 2010). While we do not yet 
understand this discrepancy, it is likely related to the 
model-predicted thermal structure. This will be a 
focus of further study in the near future. 

There are distinct differences in the predicted 
vertical distribution of dust between the two 
simulations (Figure 6). In the simulation with 
radiatively active cloud formation, the dust is more 
deeply mixed and falls off with altitude more rapidly 
than in the simulation without clouds. Preliminary 
results suggest that dust is more deeply mixed in the 
cloud case due to the radiative effects of clouds 
enhancing the strength of the mean overturning 
circulation (as discussed in Wilson et al., 2007; 
Madeleine et al., 2012). The more rapid decrease in 
dust concentration with altitude in the cloud case is 
likely the combination of two effects: clouds 
“capping” the dust by scavenging processes and due 
to increased fall velocities of particles of equal size 
with increasing altitude (e.g., see Rodin et al. 1999). 

Discussion and Conclusions: These results suggest 
that there are several feedbacks involved in 
determining the vertical distribution of aerosols in 


the Martian atmosphere. These feedbacks could 
involve interactions between the sizes of dust and 
water ice particle (and their respective fall 
velocities), dynamical responses to aerosol radiative 
effects, and possibly dust lifting processes. 

Although more work needs to be done to fully 
isolate these feedbacks, it is plausible that the key 
feedback centers on the interaction between the 
clouds and the strength and extent of the Hadley 
circulation. In this scenario, clouds that form in the 
upper branch of the Hadley cell warm the local 
atmosphere and intensify the circulation. The 
circulation then deepens, drawing more water and 
dust into the upper branch, which mixes the dust 
more deeply and allows the cloud particles to grow 
larger and interact more efficiently with the upward 
IR. The bigger cloud particles and the more deeply 
mixed dust further intensify the circulation until the 
cloud particles get big enough to fall faster, which 
limits further expansion and efficiently “caps” the 
dust at the level of the clouds. 
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